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Abstract 

We develop a resonance theory to describe the evolution of open systems with 
time-dependent dynamics. Our approach is based on piecewise constant Hamilto- 
nians: we represent the evolution on each constant bit using a recently developed 
dynamical resonance theory, and we piece them together to obtain the total evo- 
lution. The initial state corresponding to one time-interval with constant Hamil- 
tonian is the final state of the system corresponding to the interval before. This 
results in a non-markovian dynamics. We find a representation of the dynamics in 
terms of resonance energies and resonance states associated to the Hamiltonians, 
valid for all times t > and for small (but fixed) interaction strengths. The repre- 
sentation has the form of a path integral over resonances. We present applications 
to a spin-fermion system, where the energy levels of the spin may undergo rather 
arbitrary crossings in the course of time. In particular, we find the probability for 
transition between ground- and excited state at all times. 

1 Introduction and outline of main results 

We study the evolution of an open quantum system S in contact with a quantum heat 
reservoir R. The Hamiltonian of S, as well as the interaction between the two systems 
is time-dependent. Our goal is to derive the form of the reduced dynamics of S for 
all times t > and for small (but fixed) values of the coupling constant governing the 
strength of the interaction. An analysis of this kind for time-independent dynamics has 
been carried out in [19, 20]. The approach adopted in the present work is based on the 
methods developed in these references, which in turn are extensions of a recent theory of 
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quantum resonances for the analysis of large-time asymptotics of open quantum systems 
[14, 7, 11, 15, 17, 18] (see also the references in these works for further literature). 

Within the context of the present paper, the long-time asymptotics has been exam- 
ined for time-dependent dynamics in the following settings: in [12] for interaction op- 
erators having some limit as t — > oo, in [12, 2] for periodic interactions (using algebraic 
scattering theory and Floquet theory, respectively), in [9, 10] for piecewise constant 
dynamics and markovian reservoirs (repeated interaction systems) and in [1, 3] for adi- 
abatic dynamics. All these works are concerned with the approach of the system to 
an asymptotic state and with the thermodynamic properties of the latter. In contrast, 
in the present paper, we examine the dynamics of the open system for all times and 
for rather arbitrary time-dependences of the dynamics (not necessarily leading to an 
asymptotic state of the system) . 

We develop a resonance theory for Hamiltonians of the form 

H(t) = H s (t)+H R + X(t)v(t), (1) 

where H§(t) and H R are the Hamiltonians of S and R respectively, X(t) is a coupling 
constant, and v(t) is an interaction operator. We base our approach on piecewise 
constant Hamiltonians of the form (1), meaning that the Heisenberg dynamics of an 
observable A is given by 

a N {A) = e^ Hl • • • j*» HN Ae- H » HN • • • e~^ H \ (2) 

where tj > and 

H* = H 3 s + H R + \ j v j . (3) 

The dynamics (2) describes sudden changes in parameters of S and the interaction, 
and it has its own interest. A piecewise constant dynamics may also be viewed as an 
approximation of a continuous dynamics, in the appropriate limit tj — > and N — ► oo. 
We illustrate both these settings on concrete models in Section 1.2. 

The space of pure states TCs of the open system S is a finite-dimensional Hilbert 
space, and its Hamiltonian is an arbitrary self-adjoint operator on Hs- We model 
the reservoir by a spatially infinitely extended (M 3 ) gas of free Fermions in equilibrium 
at temperature T > 0. The Hamiltonian of R is given by 

H R = f \k\ 2 a*(k)a(k)d 3 k, (4) 

where the a(k) and a*(k) are fermionic annihilation and creation operators, satisfying 
the standard canonical anti-commutation relations {a(k),a*(l)} = S(k — I), see e.g. [8]. 
We understand that in (4), and for all other quantities involving R, the thermodynamic 
(infinite volume, continuous mode) limit has to be taken (this is the so called Araki- 
Wyss representation [5], see also Section A.l for further details). The equilibrium state 
of R is the quasi-free state determined by the two-point function 

^{a*{k)a{l)) = (5) 
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where (3 = 1/T, see e.g. [8]. 

The interaction operator is a sum of terms of the form 

v> = G j ® <j>{cj>), (6) 

where G J is any self-adjoint operator on Hs, and 

cf>(g>) = -L[a*(<?) + <>(<?)} (7) 

is the field operator smoothed out with a function gi £ L 2 (M 3 ,d 3 /c), called a form 
factor. Here, the smoothed-out creation and annihilation operators are defined by 



'(g) = [ g(k)a*(k)d 3 k, a{g) = [ g(k)a(k)d 3 k 



(8) 



(we take annihilation operators to be anti- linear in their arguments). Interactions of 
the form (6) induce processes of absorption and emission of quanta of R by the system 
S. 

Our approach uses a spectral deformations (generated by translation in the energy 
variable in a suitable Hilbert space). This method necessitates certain regularity of the 
form factors gK We represent g>(r,a) in spherical coordinates (r,a) € M + x S 2 and 
denote by gi its complex conjugate. 

(R) Assumption on regularity (translation analyticity) of form- factors. The maps 

R x S 2 3 (u, a) ^ a(u)\ ' ' ■ - <^ v ' J ' tc ^ n ' (9) 



gi(y/—u, a), if n < 0, 

where a(u) = l,e _/3n , extend analytically (in u) to maps from (u,a) 6{z£C : 
|Imz| < 5} x S" 2 to L 2 (M x S" 2 , difdcr), for some S > 0. 

An example of a form factor satisfying this assumption is g(k) = |/c| -1 / 2 e - ' fc ' 2 . 

The next assumption concerns the "complete splitting of resonances" . We make it 
merely for the purpose of a lighter exposition of our results. Fix j and let e G {E — E' : 
E,E' 6 spec(iig)} be an energy difference of the system S. In the resonance approach, 
the evolution of S is described by resonance energies e, which are complex in general 
and reflect the non-unitary (irreversible) character of the reduced dynamics of S. As 
the interaction between S and R is turned on, resonance energies e bifurcate out of 
each (real) energy difference e. The total multiplicity of resonance energies bifurcating 
out of a given e equals mult(e) (the multiplicity of e viewed as an eigenvalue of the 
operator Hi (g) Is — Is ® ^s ac ting on 7i.s ® TCs). We label the resonances by r = (e, s), 
where e denotes the origin of bifurcation and 1 < s < mult(e) counts distinct resonance 
energies. 

(S) Assumption on complete splitting of resonances. At all time-steps j, for Xj / 
there are mult(e) distinct resonance energies e associated to each eigenvalue 
difference e. 

One can deal equally well with degenerate resonance energies by adapting the ar- 
guments of [20] to the time-dependent case. Assumption (S) can be verified by using a 
perturbative analysis of the resonance energies (see (11) below). 
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1.1 Dynamics of S 



Let us explain our main result on the dynamics of the system S, whose precise statement 
is given in Theorem 2.1 below. We consider initial states of the form ujq = los,o ® wra, 
where u>s,o is an arbitrary state of S, and wr^ is given by (5). 1 Theorem 2.1 gives 
the following representation of the evolution of the average of an observable A of the 
system S in the initial state loq. 



u Jo {a N {A))= Y, e i Sf=i*,^(^)p rii ... ir . jv (^) + 0( m a X |A J |). (10) 



ri,...,r N 



The sum is taken over indices rj = (e,s) which label the resonance energies £ 3 {rj) 
associated to the system at step j. The p ri ,...,r N are linear functional on the algebra 
of observables 9JTg = B(Tts) (bounded operators). Both e J (rj) and p ri ,...^ N depend on 
Xj, and the remainder term depends on N as well, but is uniform in the tj > (it also 
depends on the interaction v 3 , (6)). The resonance energies have the expansion 

eP = e + X)5 3 + 0(Xj) (11) 

for small Xj. Here, the S 3 are eigenvalues of an operator A 3 (e) acting on the doubled 
space 7is <8> Hs, called the level shift operator associated to e at time-step j (see the 
definition (74)). We have ImS 3 > 0. 2 

The functionals p ri ,...,r N can be expressed as 

(A) =P(U (r u ...,r N ) A® l s ), (12) 

where P is a linear functional on the algebra 9Jts <8> Is acting on the doubled space 
<8> Hs- P depends on the initial state loq only. Ilo is a product of "transition 
amplitudes" associated to r±, . . . , rjv, 



n (ri,...,rjv) 



JV-l 

II (v S (r j ),rf +1 (r j+1 )) 

3=1 



\v l (n))(v N (r N )\. (13) 



Here, the rf(r), r/ J (r) G 7is <8> Tis are (the lowest order contributions of) resonance 
eigenvectors, defined by 



A J (eV(r) =e 3 (r)rf(r) and [A J (e)]*^(r) = ei(r)rp(r) (14) 
and normalized as 



(rf (r), rf (r)} = 1, (v 1 ( r )i 7 l' > ( r ')) = if r / r' (15) 

(recall that r = (e, s)). In (14), [A J (e)]* denotes the adjoint operator of A J (e). 

1 Our theory works as well for states which are local perturbations of such states, but we restrict 
our exposition to product initial states. 

2 This can be seen by direct calculation in concrete models, and it can also be derived from general 
considerations, see as well [18]. 
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Discussion of (10). — At each moment when the Hamiltonian changes, the system 
starts a new dynamical process with an initial condition determined by the final state 
of the previous process. Note that even if we start in an unentangled (product) state of 
S + R, already after the first bit of interaction the state will become entangled. Since 
S interacts with the same reservoir at each time-step, the dynamics is not markovian. 
The cumulative effect of the interactions is encoded in the functionals p ri ,...,r N and the 
product of increments of the dynamics e 1<l£l ^ ri - > • • • q^nz™ (r N ) ^ 

- The case of a time-independent dynamics can be recovered from (10) as follows. 
If W is close to then rf(r) is close to rf +1 (r) and the transition amplitude 

(Jf( rj ),rf + \r j+1 )) « (T?(r j ),r?(r j+1 )) (16) 

is very small unless rj = r^+i, in which case it is roughly unity (see (15)). We can view 
the sequence r±, . . . , rjv as a "path" of resonances: the system hops from resonance rj 
to resonance Tj + \ as time passes the moment t\ + . . . + tj. Thus for small differences 
Hi — iiP +1 , the main contribution to the sum in (10) comes from the constant paths 
rj = r = const, with associated propagator e lt£ ( r \ In this limit of a time-independent 
Hamiltonian, (10) reduces to the dynamics of S derived in [19]. 

- If the interaction v 3 is energy exchanging then it typically drives the total system 
S + R to its equilibrium state at a relaxation rate 1 / 7f herm . In the regime tj > Tf herm 
one then expects to find the system after each bit of constant interaction in equilibrium 
relative to the dynamics at that moment. This is an adiabatic process during which 
the state of the system follows its instantaneous equilibrium state, see also [1, 3]. 

1.2 Applications 

We consider a spin-^ particle subject to a time-dependent Hamiltonian, coupled to a 
thermal Fermi field. The space of pure states of the spin (system S) is C 2 , and the 
Hamiltonian at time-step j is given by 



A j A* 
Hi = — a z = — 
s 2 2 



1 
-1 



(17) 



where A J G M is the energy level spacing. At times j when A 3 switches its sign we say 
we have a level crossing. The interaction of S with R is given by Xv (constant in time), 
where A is a coupling constant, and (recall (7)) 



v = a x cf)(g) 



1 

1 



-L[a*(g) + a(g)]. (18) 



The parameter regime A 2 «. min,- |A J | describes well separated resonances, while if A 2 
is of the order of the , or if A 2 » min,- 1 A 1 \ , then we have "overlapping resonances" . 

1.2.1 Regime of overlapping resonances 

Our goal is to analyze the dynamics of the system for small but independent values 
of A and A J . In Theorem 3.1 and Proposition 3.2 we give the explicit form of the 
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resonance data e 3 (r),rf {r),rj 3 (r) for this model, as well as the transition amplitudes 
(rp(r), r/ J+1 (r / )) (recall (13)). As an illustration of these results, we present here the 
case of a single sudden level crossing and the limit of continuous time-dependent dy- 
namics. 

Single sudden level crossing. Consider a single sudden level crossing at time 
t c , whose evolution is generated by (18) and 

H f s = ^a z , 0<t<t c , H f s = -^a z , t c <t, (19) 

where Ai, A2 > 0. Denote by p ge (t) the probability that the system S is at time t in 
the excited state of while at time zero it started off in the ground state of -ffg =0 - 
(By ground state we mean the state with lowest energy.) We show in Section 3.1.1 that 
for independently small values of A and A max = max{Ai, A2}, we have 

If l-e ite if0<t<t c , 1 a \ , onA 

Pge (t) = - 2 [ 1 + eite ift " tc) )+0(|A| + A max ), (20) 

where 

e = i 7r A 2 7o + 0(A 2 (|A| + A max )) (21) 

with 

70 = lim f / da\g(^v)\ 2 - (22) 

r^0+ 2 J S 2 

It is assumed here that < 70 < 00, which amounts to an infra-red (\k\ ~ 0) behaviour 
g(\k\,a) ~ |A;|~ 1//2 in three space dimensions (spherical coordinates; see also Assumption 
(R) in Section 1). Formula (20) shows in particular that at t c the probability jumps 
up by an amount 

5 = e iu - £ + 0(|A| + A max ) = e -7oA 2 t c [i+o(|A|+A max )] + G( | A | + Amax) (23) 

This is in part explained by the fact that the excited state itself jumps at t c from |+) 
to |— ). As t — > 00, Pge(t) approaches 1/2 + 0(|A| + A max ), which is the probability of 
finding S in the excited state when the total system S + R is in equilibrium, provided 
A and A max are small. This is the correct value of this asymptotic probability, since 
the system exhibits return to equilibrium. 

In a time-dependent setting where energy levels of a quantum system are brought 
close together (but do not cross) in the course of time, say due to some external forc- 
ing, a transition from one energy state to another is called a Landau-Zener transition. 
Landau-Zener theory is important in physics and chemistry, see e.g. [13, 21]. The influ- 
ence of dissipation on Landau-Zener transitions has been studied for different systems 
and transition probabilities similar to (20) have been calculated [4, 23]. We mention 
that our method allows for crossing of the energy levels in the course of time. 

Continuum limit. We investigate the continuum limit of the model (17) by 
setting tj = jjj, j = 1, . . . , N, with fixed t and letting iV — > 00. Our goal is to find the 
limit of the dominant contribution to the dynamics, given by the sum in (10). We do 
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not control the remainder term in (10) uniformly in N in this paper (work on this is 
in progress). 

We take A J = A(jjjj), where A(t) is a continuously differentiable function of t > 
(of course, one can deal with less regular A if desired), and we define r- 7 and a by 

= —, a = ^A 2 70 , (24) 

where 70 is given in (22) above. For simplicity of the exposition, we will assume in 
what follows that < r max < 1, where r max := A ™ ax and A max = sup t>0 |A(t)|. This 
regime is interesting since it accommodates the situation of level crossings (A(t) = 0) 
while A is fixed. 

We have seen above (see (16)) that if A J+1 — A J is small, then the transition 
amplitude (rp(r), rf +1 (r')) associated to a jump (r / r') in the resonance path at j is 
small. We show in the proof of Theorem 3.4 that this amplitude is at most of the size 

W:=sup|r'(t)|. (25) 
t>o 

(The ' here means derivative.) The sum over all paths in (10) can be written as a sum 
over all paths with k jumps, k = 0, . . . , N — 1. This sum becomes an infinite series in 
the continuous time limit, and the summand associated to a path with k jumps is of 
the order of (r max ) fc . For small r max , one can thus (rigorously) approximate the series 
by the first few terms. In Theorem 3.4 we show that the continuous time limit of the 
sum in (10) is given by 

4 

J2 ei/o ^)d^ (0 , t, r) (Vo, B V (0, r)> {rj(t, r),A^ ci ) (26) 

r=l 

Jo i + y+(s) 2 

X <V>0,£??(0,3))<^(MMY'ref> 

+ /** e 1 ti'V, 4)d s ' + i/ s * e(s ' ) 3)d s ' w( p t8i 4) v -W+W w (8,t,3)d3 
Jo l+y^{sY 

X (^ ,^(0,4))(^,3),^Vref) 

+ o([t^]V^). 

Here, e(t, r) and rj(t, r), rj(t, r) are the resonance energies and resonance vectors at time 
t > (with r = 1, ... ,4; see Section 3.1.2 for their explicit form). The functions w 
and y± are associated to the transition amplitudes, ipo is the (Gelfand-Naimark-Segal) 
vector representative of the initial state of S (represented in the Hilbert space Tis'&'Hs), 
V> re f is the vector representing the trace state of S, and B G Is <8> 9Ks is the unique 
operator satisfying Vo = ^Yw- We refer the reader to Section 3.1.2 for the explicit 
expressions of all quantities involved in (26). 

The sum in (26) is the contribution coming from the constant paths, while the 
two integrals come from paths having a single jump (taking place at the integration 
variable s). Both these integrals are of the order of r max . Naturally, the Riemann sum 
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in the propagator of (10) has now become an integral, and the product of the transition 
coefficients has turned into the factors w (associated to products without jumps) and 
the fractions involving the y± (associated with single jumps). The remainder term 
contains the contributions of all paths with more than one jump. In fact, it is not very 
hard to find the continuous time limit of all paths with arbitrarily many jumps. We do 
not present the corresponding formulas in this work since they are rather cumbersome. 



1.2.2 Regime of isolated resonances 

The regime of isolated resonances can be obtained as a limit (t j » 1) of the case 
of overlapping resonances, but not the other way around. More specifically, the ex- 
pressions for the resonance energies and resonance vectors obtained for overlapping 
resonances (see Theorem 3.1) are still valid once equations (22) and (24) are replaced 
by 



7(A'') = 
a(A j ) = 



da 



S 2 



AJ\,a 



A 2 7 (A^'). 



(27) 
(28) 



We present the explicit form of the resonance data in Theorem 3.5. 

Single sudden level crossing. With the same set-up as before, but now assuming 
that A 1 , A 2 > are fixed independently of A, 



Pge(<) 



1 _ g-^AtiCA 1 ) 

e^ 1 + 1 



+ 0(|A|) 



(29) 



for t < t c and 



Pge(*) 



1 



-7rA 2 t c 7(A 1 ) 



+ 



-7rA 2 (t-i c ) 7 (A 2 ) _ 



-/3(A 1 +A 2 ) 



e"^ 2 + 1 



+ 0(|A|) (30) 



for t > t c . Therefore, there is a jump up at t c equal to e ^-^(A 1 ) _|_ 0(|A|), in perfect 
analogy to the situation before. 
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2 Dynamical resonance theory 

We give a precise definition of the model in Section 2.1. In Section 2.2 we prove formula 
(10). The main result is Theorem 2.1. 
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2.1 Description of the model 

The Hilbert space of states of the system S + R is given by 

H = H s ®Hs®F, (31) 
where Tis = C d is the Hilbert space of pure states of S, and 

T = F(L 2 (R x S 2 ,du x dE)) (32) 



is the fermionic Fock space over the one-particle space L 2 (R 3 x S 2 , du x dE). 7i is the 
(Gelfand-Naimark-Segal) representation Hilbert space associated to the reference state 

^rcf = ^S,ref ® Wr,/3- (33) 

Here, ^s.ref is the trace state on the C* algebra of observables 2ls = B(Hs) (bounded 
operators on TCs), and wr^ is the equilibrium state of the infinitely extended free Fermi 
gas (see also (5)) on the C* algebra 21r generated by the creation and annihilation 
operators {a(g), a* (g) : g G L 2 (IR 3 , d 3 k)} (called the Canonical Anticommutation 
Relation (CAR) algebra). We refer the reader to Appendix A. 2 for more information 
on this representation of the CAR algebra. 

The Hilbert space (31) supports in particular all local modifications to of uj xc {. Such 
an lu has the form 

uj(A) = <^,7r(AhU, (34) 

for all A 6 21 = 2ls ® 21r, for some ij) w G H, where 7r : 9Jt — > B(7i) is the representation 
map. 

The dynamics of the system generated by iiP, (3), is represented on TL by a Liouville 
operator U : 

n (jtHi Ae -itWj = jtU^ A)e -itU (35) 

for all A G 21 and all < £ R. Consequently, (2) is represented on 7i as 

e itlLl • • • e i ' iviiV 7r(y4)e- i * JvLiV • • • e~ itlLl . (36) 

(Of course, we understand that the thermodynamic limit has been performed.) The 
Liouville operators have the form (see Appendix A. 2) 



where 
Here, 

acts on Hs <8> and 



V = L{ + XjV j , (37) 

4 = 4 + L R . (38) 
4 = H J S ® 11 - 11 (g> ^ (39) 



L R = dr(n) (40) 

is the second quantization of multiplication by a £ 1 acting on T. The interaction 
operator 

V j = n(v j ) (41) 
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belongs to the C* algebra vr(2l) . 

It is useful and standard to consider the weak closure of 21, 

Wl= (2l s ®2l R )" = £(Hs)® l S ®2l R . (42) 
SDt is a von Neumann algebra acting on H. We introduce the reference state 

V'ref = tpS ® tpR, (43) 

where ips £ <8> represents the trace state of S, and where V>R G ^ is the vacuum 
vector of J 7 , representing the equilibrium state of R. The vector ip Te f is cyclic and 
separating for the von Neumann algebra 9Jt, and we denote by J = Js (g> Jr, A = 
Is <8> Is <8> the modular conjugation and the modular operator of the pair (9JT, ^ref ) 
(see also [8, 15]). It follows from the form of the interaction, (6) and Assumption (R) 
in Section 1 that A^WA -1 / 2 G 971 for all j. 

2.2 Proof of (10) 

Our main result on the piecewise constant dynamics is the following. 

Theorem 2.1 (Dominant paths) There is a constant c > s.t. if max.,- |Aj| < c, 
then we have the following. Let A £ 9JTs be any observable of S, and let tpo be any 
initial state of S + R, given by Bip Te { some B G 9JTg. Then 

(^o,e lhLl • • • jtNL N Ae -it N L N . . . Q-ihL 1 ^ ( 44 ) 
= J2 e 1 ^ t ^^(^0,BU (r 1 ,...,r N )A^ Te{ )+O(ma X \X j \), 

3 

n,...,r N 

where the e 3 (r) are the resonance eigenvalues (see (11) and also Propositon 2.7). The 
error term in (44) is uniform in the tj > 0. Let rf(r), rp{r) be the resonance eigen- 
vectors (see (14) and also (73)). Then 

n (ri,...,rjv) 

= IvHn)) (vHnWiri)) ■ ■ ■ (v N -Hr N ^), V N (r N )) (v N (r N )\. (45) 

Remarks. 1. We think that a more detailed analysis of the remainder term in (44) 
would yield an estimate 0(maxj \ Xj\) uniformly in N, but we do not prove this here. 
2. Theorem 2.1 implies formula (10). 

The remaining part of this section is devoted to the proof of Theorem 2.1. We build 
up the proof in several steps. 

2.2.1 Passage from the operators LP to the operators K 3 

Let (V'y be any operator belonging to the commutant 9JT', 3 and set K 3 = L? ' + Xj(V') 3 
(with domain V(K 3 ) = V{L 3 ) = V(L,{)). We define the operator e itKJ , t £ R, via the 

3 The commutant of SOT, (35), is defined as W = {A £ B(H) : AB = BA \/B € SOT}. 
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operator-norm convergent Dyson series 

e itKj = Va* /V--- r~ 1 ds k (v'y(s k )---(v'y(s 1 y tLj , (46) 

fc>0 - 70 ^ 

where (V') j (s) = Q lsL ' {V'y e~' lsLj . 

Since JA 1 / 2 Aip rc{ = A*ip rci for all A G 9Jt, and since V"-? = (V- 7 )*, we see that upon 
choosing 

(V'y = -JA^VijA 1 / 2 , 

we have iP^ re f = and hence 

e ltKJ ip rei = ip rei . 

Note that JA 1 / 2 = A" 1 / 2 J, and that J9JIJ = SDT', so that (V') j G 97? indeed. 
Proposition 2.2 M^e /taue 

e it^ Ae -i^ = jtLi Ae -itL^ (4?) 

/or all t £ 1, A £ 9Jt. 

Proof. It is easy to verify directly that for <f> G T>(Lq), we have ^e ltK3 4> = e ltKJ (f>. 
We write (46) as e 1 *^ = S'(t)e ltL3 = e ltLJ R'{t), where 5"(t) is given by the series on 
the r.h.s. of (46), and R'(t) = e~ itLj S(t)e itL3 . Both S'{t) and R'(t) belong to the 
commutant Wl'. Consequently, we have for all A G 97T 

e^'Ae - "^ = S'(t)d tLj Ae- itDi R!(-t) = e ltL3 Ae~ itL3 S'(t)R'(-t). 

Furthermore, S'(t)R'(-t) = e itR3 e~ itL3 e itL3 e~ itK3 = 1, and thus we obtain (47). □ 

2.2.2 Resolvent representation of propagators 
Proposition 2.3 Let A G Wl and ip G H. We have fort>0 

(V, e HKi A^) = =± jf e i2t (V, - z)- l A^ ci ) dz, (48) 

ifj> C\Xj\ for some C > 0. 

Proof. The function t ^ (ip, e ltL3 Ae~ ltLJ ip re f^ = (ip, e ltR3 Atp TC ^ is bounded and 
continuous in t G R. It follows that 

r7+ioo 



W tKi A^J) = ^-j [e 2t y o e— (v,e i ^^ ref )d S Jdz 

i />7+ioo 

= — / e zt (V, i(Ki + iz)- 1 ^ > dz 

/ e i2i (V, (K J - z)- 1 ^) dz. (49) 



-1 

2vri 



t— 17 

In the first step, we use the Laplace inversion theorem (see e.g. [22], Chapter II, 
Theorem 9.2) and in the second step we integrate the propagator to obtain the resolvent. 
□ 
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2.2.3 Analytic continuation of matrix elements 

For 9 G R we define the unitary group Ug = e 9dr ^ 9u ^ on J 7 , (translation in the radial 
variable u, see (32)), and we set 

Ll(0) = UeLlUe 1 (50) 

An easy calculation gives L J (6) = L J S + Lr + ON, where N = dT(l) is the number 
operator on T ' . Accordingly, we define in the same way 

k\q) = v{e) + \ j P{e), (51) 

p(0) = V j {0) + {V') j {0). (52) 

Let > be fixed and define the strip 

S 0Q = {9GC : \lm9\ < }. (53) 

In accordance with analytic spectral deformation theory, we assume the following an- 
alyticity condition. 

(CI) i— > P (6) has an analytic continuation as a map from Sq to the bounded oper- 
ators on H, and sup 6eSe \\P (0)\\ = C < oo. 

Proposition 2.4 If the form factors gi satisfy Assumption (R) in Section 1, then 
Condition (CI) is satisfied. 

The proof of this proposition is not hard (one examines directly the explicit expres- 
sion for P{9), see also Appendix A. 2). The following is the key technical result of the 
analytic deformation method. 

Proposition 2.5 Take z with Imz < — C|Aj|, where C is the constant in Condition 
(CI) above. The map 9 ^ [K^{9) — z)^ 1 has an extension from 9 G M to < Im# < #o- 
This extension (denoted by the same symbol) is an analytic map from {9 G C : < 
Im(9 < #o} to the bounded operators ofTL, and it is continuous in the operator norm as 
lm9 i 0, at all 0^0. 

Proof of Proposition 2.5. We fix the index j and omit it from the notation. Let 9 G 
R. On£>(L )n£>(A0 we have L {9) = L +9N and so, by (CI), K{9) = L +9N+XI{9) 
has analytic continuation to 9 G Sg (as a closed operator on V(Lq) n V(N)). The 
spectrum of the normal operator Lq + consists of real eigenvalues e G spec(Ls) and 
of horizontal lines {n9 + R : n G N} of continuous spectrum. (Note that Lq + 9N is 
the sum of two commuting self-adjoint operators.) For Im# / the eigenvalues e are 
isolated. It follows from standard perturbation theory that the spectrum of K(9) lies 
within a distance of |A| ||/(#)|| from that of Lq + ON. 

For < Im# < 0q, all Imz < — C\X\ (with C as in condition (CI)) belong to the 
resolvent set of Lq + ON, as well as to the resolvent set of K(9). For such z we express 
the resolvent using the norm-convergent Neumann series 

(K(9) - zy 1 = (L + 9N- z)' 1 ^(-X) n [I(0) {Lq + 9N- z)" 1 ]"- 

n>0 

It follows that for all z with Imz < — C\X\, 
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1*. N(K(9) - z)- 1 is bounded for all nonzero 9 with < Im0 < 9q, 

2*. (K(9) - z)- 1 is bounded uniformly in 9 s.t. < lm9 < 9 , 

3*. Ran(if (0) - z)- 1 C V(L ) n V{N) for all nonzero 9 s.t. < Im0 < 9 . 

Let 9 be s.t. < Im9 < O and take A9 to be small so that < lm(0 + A9) < 9 . The 
resolvent identity gives 



(K(9 + A9)-z)- 1 -(K(9) 



z 



-i 



{K{9 + A9) - z)- 1 [A9N + XI{9 + A0) - XI(9)} (K(9) - z)' 1 . (54) 

Relation (54) (together with the above points 1* and 2*.) shows that 9 i— > (K(9) — z)~ l 
is continuous on {9 G C : < Im0 < 9q} in the topology of bounded operators. 
Furthermore, 

(K{9 + A9) - z)- 1 - (K(9) - z)- 1 
A9 

= (K(9 + A9) - z)- 1 [N + XX(9, A9)} (K{9) - z)' 1 , (55) 

where limAe^o X{0, A9) = del{9). Combining (54), (55) and points 1*, 2* above, we 
see that for all Imz < — C|A| and 9 s.t. < Im0 < 9q, 

d e (K{9) - z)- 1 = (K(9) - z)' 1 [N + Xd e I(9)] (K{9) - z)~\ 

the r.h.s. being a bounded operator. 

We now show that (K(9) — z)^ 1 is continuous as0 = x + iy^xG M\{0}, y > 0. 
The resolvent identity gives 

(K(x + iy)- z)- 1 = 

{K{x) - z)- 1 + {K(x + iy) - z)- 1 [K{x) - K(x + iy)] (K(x) - z)~\ (56) 

We rewrite (56) in the form 

(K(x + iy) - z)~ l [t - W{y)\ = (K(x) - z)~\ 

where 

W(y) := [K(x) - K(x + iy)){K{x) - z)' 1 — ► 
in operator norm, as y —>■ (here it is important that x^O). It follows that 



lim (K(x + iy) - z)' 1 = lim (K(x) - z)' 1 ^ - W^y)}- 1 = (K(a 



This completes the proof of Proposition 2.5. □ 
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2.2.4 Separating pole contributions 

For A G SDts and < j < N - 1 we set 

Aj = e ^+ lL3+1 • • • jtNL N Ae -it N L N e -it J+1 v+i ^ ( 57 ) 

and we define An = A. Let -B be the operator in the commutant 9Jtg satisfying ipQ = 
.BVref) and denote by P^ R = |V'r)(V ; r| the orthogonal projection onto TCs <8> Ti-s <8> C?/>r- 
If all Aj = 0, then the dynamics of S and R decouple, and we have Aj G 9#s ® 1r> 
and thus Ajtp rc f G RanP^, R . The following result follows from an easy perturbation 
expansion. 

Lemma 2.6 Sei i^, R = 11 — P^, R . PFe Ziawe 

lli^, Aj^refll < Cjmax|Aj| ||,4||. 

j 

Using this result, we arrive at the representation 
, ^ ,e ihLl A 1 e- ihLl ^ ) = ^B^A^) = (^0, Be ihRl P^A^J) + R u (58) 

where ||i?i|| < Cmaxj \Xj\ (with a constant C depending on N and ||^4||). The scalar 
product term on the right side of (58) can now be treated as 

[ih,B^ Kl P^A^J) = (ij ,Be ihKl P^e it2K2 A 2 ^ 

^o, B^ Kl P^ K2 P^A 2 ^ rei ) + R 2 , 

with ||i?2|| < Cmaxj \Xj\. We iterate this procedure and arrive at 
, . 0> ^ . . . jt N L» Ae -it N L» . . . e-i^ 1 ^ 

> , BP^ Kl P^ ■ ■ ■ P^J tNKN P^A^ rc[ ) + i?, (59) 

where the remainder term 

R = ^Rj (60) 
j'=i 

satisfies ||i?|| < Cmaxj |Aj|, with C depending on N. 
Using the resolvent representation (49), we obtain 

Pi>J tjKj Pi>n = ^ / e^Ji^ - z)- 1 ^^ dz. (61) 

We now perform spectral deformation in the integrand on the right side of (61). For 
SeKwe have P^Ug = P^ R = UgP^ R and it follows that 

P^{K j - z)-^ = P^{K\9) - zr'iV (62) 
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Proposition 2.5 shows that (62) has an extension to values of 9 in < Im9 < 9q, and 
that this extension is analytic in the open strip < lm9 < 9q, and continuous on 
R\{0}, provided Imz < — Cmaxj \Xj\. However, since (62) is constant for 9 on the real 
axis, it must actually be constant on the entire region {0 < lm9 < #o}- 4 This shows 
that (62) holds for Imz < —C max.,- |Aj| and < lm9 < 9 . We thus have 

P^ K3 P^ = =± f &*P^(Ki(6) - z)- l P^ dz, (63) 

for all 9 with < Im# < 9q , and where 7 > C maxj | Xj \ . We analyze the integral on 
the r.h.s. of (63) in more detail. The following is a standard result [19, 20, 14]. 

Proposition 2.7 Fix 9' with < 9' < 0q. There is a constant cq > s.t. if |A| < cq/(3 
then the spectrum of K^{9), lying in the complex half-plane {z G C : Imz < 9'/2}, 
is independent of 9 in the region 9' < lm9 < 9q. It consists of the distinct isolated 
eigenvalues 

•{V(e, s) : e G spec(Lg), s = 1, . . . , f(e)}, 

where 1 < u(e) < mult(e) counts the splitting of the eigenvalue e into distinct reso- 
nances. Moreover, we have e J (e, s) — > e as A — > 0, for all s, and Ime(e, s) > 0. The 
continuous spectrum of K^(9) lies in the region {Imz > 39' /4}. 

We now "shift" the path of integration R — of the integral in (63) to the path 
R + 3i6*'/4 in the upper half plane. Thereby we pick up contributions (residues) coming 
from the poles of the integrand. Let C J (e,s) be a small circle around e J (e,s), not 
enclosing or touching any other spectrum of K 3 (9). Define the generally non-orthogonal 
Riesz spectral projections 

Q 3 (e,s) = ^f {K^9)-z)- l Az, (64) 

and the operator 

Qj( 00 )=Qj( 00: t j ,9) = — [ e^{K\6)- Zj^dzj. (65) 
For any vectors tp, 4> G H we have by standard contour deformation of complex integrals 
-± / ^{iP,{K\9)-z 3 y^)dz 3 (66) 

„(e) 

= E E ei ^ (e ' s) + <^> & (°°)^> • 

eSspec(L^) s=1 

The operator Q J (00) reduces to e ltjLo ^P^ R for Aj = 0, and one can show the following 
result. 



4 Apply the Schwarz reflection principle to the analytic function F{9) — F(0), where F{&) equals 
(62). 
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Proposition 2.8 ([19]) We have \\P^ R Q j (oo)P^J < CA]e~ 3t ^'/4_ 
Combining relations (59), (63), (66) and Proposition 2.8, we obtain 

e^ 1 . . . jt N L» Ae -u N L» . . . e-i^ 1 ^ (67) 

= e i ^ t ^^)(B>o,n(r 1 ,...,r^ r ef) + i? / , 

ri,...,r N 

where the (multi-)indices rj are summed over {(e, s) : e £ spec(Lg), s = 1, . . . , ^(e)}, 
with e^(rj) = e J '(e, s) if r,j = (e, s) and where we introduce 

U(r 1 ,...,r N ) = Q 1 (r 1 )---Q N (r N ), (68) 

with Q J (rj) given by (64). The remainder term R' in (67) satisfies ||i?|| < C maxj |Aj| 
(with C depending on N). 

As explained in the introduction, we assume that 

(S) Each projection Q^(e, s) has rank one (for Xj > 0). 

This assumption means that all resonance energies £ J (e, s) are simple, and it is valid in 
all our applications. One may modify the results of [20], where degenerate resonance 
energies are treated for time-independent Hamiltonians, to eliminate Condition (S). 
Having rank one, the projections are given by 

Q i (r j ) = \x i (r j ))(x j (.r j )\, rj = (e,s), (69) 

where 

K j (e) X j (rj) = ei( rj )x j (rj), (70) 

= [^)]*X%), (71) 

and the resonance eigenvectors are normalized as 

(x j (r j ),X j (r j )) = l- (72) 

Using perturbation theory (e.g. the Feshbach technique, [6, 19, 20]), one sees that 
the resonance eigenvectors have the expansion 

X j fa ) = rf fa ) <8> ^r + 0(Xj) and x j fa ) = 7p fa )®^R + 0(\j), (73) 

where Vr is the vacuum vector of P, and where the vectors J? 7 fa) and fpfa), for 
Tj = (e, s), belong to the eigenspace of L J S associated to the eigenvalue e. Let Pi be 
the orthogonal spectral projection of Lg associated to the eigenvalue e. We define the 
level shift operator A- 7 (e) by 

AJ ( e ) = pop pi ( e ) (p> Q - e + iO)" 1 P j {e)P P 3 e , (74) 
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where Pi = 1 — Pi and where Lq = PlL J Pl \ RanPi The vectors if(e,s) and 
rp(e,s) are eigenvectors of A 3 (e) and its adjoint [A J (e)]*, 



A J (e)rf(e,s)=5 J (e,s)rf(e,s) and [A J (e)]*rf (e, s) = 5 3 (e, s)r/ J (e, s), (75) 

satisfying the normalization relation 

(rf(e, S ),n i (e,s)) = l. (76) 

The resonance energies have the expansion 

£ j(e, s ) = e- X]P{e, s) + O(Xj). (77) 

The proof of Theorem 2.1 is now complete by combining expansion (73) with (67) and 
(68). 

3 Applications: details and proofs 

The setting of the applications is given in Section 1.2. 
3.1 Regime of overlapping resonances 

The explicit (perturbative) form of the resonance data for the system is given in the 
following theorem. Recall that t 3 and a are defined in (24). 

Theorem 3.1 (Resonances of K) There is a constant C s.t. i/|A|+maxj \A 3 \ < C, 
then the resonances of K 3 in the region {z G C : Imz < 9q} ar & given by 

e j {l) = (78) 
eP{2) = 2kr + 0(A 2 (|A| + A max )) (79) 

eJ(3) = i<7 + o^j (ri)2 - 1 + 0(A 2 (|A| + A max )), (80) 

eJ(4) = id - a^(ri)2 _ \ + (\ 2 (\X\ + A max )), (81) 

where the square root always means the principal branch, with branch-cut on the negative 
real axis (the argument function takes values in (— ir, tt]). The resonance eigenvectors 
X J (r) and x J ( r ); r = 1) • • • ,4, are given by (73), with 

7/(1) = rP(l) = ^=[<p ++ + <p-] (82) 

rf{2) = r7(2) = -L^ ++ -v9__] (83) 

^'(3) = v+ .+y 3 +t p.+ , ?? {3) = a? + (<p+- (84) 
r] j {4) = y + _+y>_ip- + , rp'(4) = a>_(ip + - + y j _ip-+), (85) 



where 



y± 



i = -i^±i^/(ri)2-l, «i = [l + ( 1 4) 2 ]- 1 . (86) 
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It follows from Theorem 3.1 that all resonances are non-degenerate provided (t-' ) 2 7^ 
1 for all j, a condition we assume to hold in this section. Hence condition (S) is satisfied 
(see introduction as well as Section 2.2.4). We define the transition coefficients 

T3(r,r') = (^(r), V i +1 (r')), r, r' = 1, . . . , 4, j = l,2,..., (87) 

in terms of which the dominant path Ho(ri, . . . , rjv), (45), can be written as 

N-l 



U Q (n, ...,r N ) 



^(n))^^)!- (88) 



Before proving Theorem 3.1, we mention that expressions (82)- (86) yield the following 
result. 

Proposition 3.2 (Transition coefficients) We have 

T j (l,l) = T J (2,2) = 1, 

2* (3, 3) = l+^yi[yi +1 -yi], 

T'(4,4) = l + Zyi[yi +1 -yi], 

T'(3,4) = 



^4// oi/ier transition coefficients vanish. (Note that if A^ +1 = A J then T J (r, r') = 5, 
(Kronecker symbol)). 

In the regime of separated resonances, where r m i n := min^ |r J '| » 1, we have 



r.r 



T'(3,3), T'(4,4) = 1+ — T J ^^ + 0(1/^), 



T»(3,4), T>(4,3) = — ± J— ^ + 0(l/r 2 in ). 

In £/ie regime of overlapping resonances, where r max := maxj \t j \ « 1, we /ictue 

T'(3,3), T'(4,4) = l±i + 0(r r 2 nax ), 



2 

T'(3,4), T'(4,3) = ±i + 0(r 2 ax ). 



Proof of Theorem 3.1. Throughout the proof, we consider j fixed and do not 
display it. The unperturbed Liouville operator Lq = Lr has a four-fold degenerate 
eigenvalue at the origin, and absolutely continuous spectrum filling the entire real axis. 
We have K = Lq + I, where 

I=^L S + X[V-V], (89) 
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with Ls = cr z <8> 1 — 1 <8> & z , V = a x <8> Is <8> Vpig) and V' = Is <8> cr x ® <pp(g)- The 
spectrally deformed operator, for G 5e , is given by 

tf(0) = L o + #iV + /(0), (90) 

where N = dT(l) is the number operator in J 7 , and where 

m = ^L s + \[V(9) + V'(9)}. (91) 

We consider 9 = W , with 9' > 0. Let P be the projection onto TCs <8> ~Hs <8> C^r, and 
let P = 11 — P. We denote by T the restriction of an operator T to the range of P . 
The estimate || P (L R + i0lV - z)- l \\ = [dist(z, ifi'N* + K)]" 1 implies 

llP^R + i^'iV-z)- 1 !! = ^, forlmz<#72. (92) 

The Neumann series 

P(K(i9') -z)- l P 

= PiLn + W'N-zr'P^i-^rUi^iLR + ieN-z)- 1 ^ (93) 

n>0 

converges for all Imz < 9'/2, provided that ||7(i0')||^ < 1. The latter condition is 
satisfied for |A| + |A| < C(6*o), see the assumptions in Theorem 3.1. We consider 
spectral points z with Imz < 6'/2. The Feshbach map method [19, 20, 6] tells us that 
such a z belongs to the spectrum of K{\9') if and only if it belongs to the spectrum of 
the Feshbach map applied to K(i9'), 

F PtZ {K(\9')) = P [K{\9') - I{\9')P(K(\9') - z^PIiiO')} P. (94) 

Using (91) and (93) we obtain 

F PiZ (K(W')) = P^L S (95) 

-\ 2 P[V{i9') - V'(i9')] P (L R + W'N - zY l P [V(i9') - V'(W')]P 
+0(A 2 (|A| + |A|)). 

Furthermore, the estimate (L R + i9'N - z)~ l P = (L R + i9'N)- l P +0(\z\/(9') 2 ) leads 
to (fixed 9') 

F P , Z (K(W)) = A + 0(A 2 (|A| + |A| + |s|)) , (96) 



where 



A = jL s - A 2 A R (97) 



is the level shift operator. Here, 

Ar = P[V - V] P (I R + iO)- 1 P[V- V']P (98) 
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We understand (96), (97) and (98) as operators acting on RanP = C 2 <8> C 2 . In (98) 
we have eliminated the spectral deformation parameter W by analyticity in a standard 
fashion, replacing it by the (operator norm) limit of (Lr + ie) _1 P as e J. 0. 

Our next task is to calculate the eigenvalues and eigenvectors of the level shift 
operator A. Using the explicit form of V and V we obtain the following result by a 
standard and straightforward calculation (see also [16, 9, 19, 20], for instance). 

Lemma 3.3 We have A 2 Ar = — ia + ia (a x a x ), where a = ^A 2 7o, with 70 given in 
(22), and where a x is the Pauli matrix (c.f. (18)). 

In the ordered orthonormal basis {ip ++ ,ip^ ,ip ^,ip } of C 2 ® C 2 , the level shift 

operator (97) is represented by the matrix 












— ia 





A 


— ia 








— icr 


-A 





— ia 












A-1T+ ; t \ : 

—ia 

It is now a simple matter to verify that 

spec(A) = {0, 2icr, ia[l T V(AAr) 2 - 1]}, 

with corresponding eigenvectors r?(l), • • • ,^(4) given by (82)-(85) (with j fixed). 

It follows from (96) and the isospectrality of the Feshbach map (mentioned before 
(94)) that the resonance eigenvalues e(r) are given by (78)-(80). 

One equally easily finds the eigenvectors rj(l), . . . , 77(4) of the adjoint of (99). 

This completes the proof of Theorem 3.1. □ 



3.1.1 Single sudden level crossing 

The Hamiltonian is given by (19). The following is a proof of the expression (20) for 
the transition probability p ge (t). 

The probability p gc (t), for t > t c , is given by 

Pge(t) = (ip ,j tcLl e i{t -^ L \A Is ® HrK^-'^V^Vo) , (100) 

where -00 = f ® ipR, and where A = |<^_)(<^_|. (Note that (p- is the excited state of 

Hg.). We use Theorem 2.1 with B = \/2 Is <8) \(f-)(ip-\ (so that Bijis = V > where 

ipS is the trace state of S). The sum in (44) has only two terms (since we have only 
one jump in the Hamiltonian), and we obtain 

Pge(t) (101) 

= e i *<= £l ^) +i (*-*c)- 2 (^)r(r 1 ,r 2 ) <^__,r ? 1 (r 1 )> (rf(r 2 ), <^__> + 0(|A| + A max ), 

ri,r-2 

where the remaider here contains an 0(A max ) term since we carry out perturbation 
theory in the coupling constant A and the energy spacing A simultaneously (overlapping 
resonance regime). Theorem 3.1 and 3.2 then imply that 

p gc (t) = \ + ^(2 )+i (t-t c )eH2) + om + Amax); 
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which shows (20) for t > t c . The proof for < t < t c goes along the same lines (and is 
actually easier, since there is no jump in the Hamiltonian in this case). □ 

3.1.2 Continuum limit and slow variation expansion (r'(i) small) 

We investigate in this section the continuum limit, as explained in Section 1.2. 

All quantities in Theorem 3.1 and Proposition 3.2 can be viewed as depending on 
continuous time t by the substitution j i-> jjj. For instance £ J (3) = £(jj?, 3), where 



e(t, 3)=ia + iay/1 - r(t) 2 + 0(A 2 (|A| + A max )), (102) 
with t 3 =r{jjj) and y 3 ± = y±{jj r ) with 



y±(t) = -ir(t)T V^~r(t) 2 . 

Similarly, the continuous time resonance vectors are denoted by r](t,r), rj(t,r); for 
instance, 

»7(*> 3) = <p+- +y+(t)(p-+, rj(t,3) = a + (t)(<p+- +y + (t)tp-+), (103) 

and so on. 

We define for s, t > the quantities w(s, t, 1) = w(s, t, 2) = 1, and 



It is quite clear (see the proof of Theorem 3.4 below for details) that the transition 
coefficients T J (r, r') associated to a jump (r / r 1 ) are of the size 

r; ax :=sup|T'(t)|. (105) 

Now the sum over all paths in (44) can be written as a sum over all paths with k 
jumps, k = 0, . . . , N — 1. This sum becomes an infinite series in the continuous time 
limit. Each path with k jumps is of the order of (T4, ax ) fc , and so for small t^ x , one 
can approximate the series by the first few terms. 

Theorem 3.4 (Continuous time limit) Let r](t, r) and rj(t, r) be the resonance vec- 
tors at time t > 0, where r = 1,...,4 (see (103)). The continuous time limit of 

the main term of (44), J2 ri r N e 1 ^ 1 * j£ ' (rj) (B*tp , U (n, . . . , r N )Atp rei ), is given by 
(26). 



Proof of Theorem 3.4- By the mean value theorem we have 



y± +1 -v± = x 



-ir'(t 1 )±r'(t 2 ) T(t2) 



aA-t^) 2 
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where ti,t 2 G (jt/N, (j + l)t/N). Thus the transition coefficients satisfy 



T*(r,r) = l + 0(r max t/iV), Vj,r (106) 

T ' (3 ' 4) = 7f 1 + yXvN) [y ' +{jt/N) + ° (rLx t/N)] ' (1 ° 7) 

and T J (4, 3) is given by the r.h.s. of (107) with y + replaced by y_. The remainders 
0(T^ iax i/iV") are uniform in j and r. Relations (106), (107) show that for "slow varia- 
tions" of r(t) the transition coefficient associated to a jump 3 <-> 4 is small, proportional 
to T'(t)t/N, while no-jump transitions have weight one. 

We write the sum over all paths in (44) as a sum over all paths with exactly k jumps, 
where fe = 0,l,...,iV — 1. (A jump happens if one value of rj changes to a different value 
of rj+i). There are exactly four paths without any jumps, corresponding to IIo(r, . . . , r), 
r = 1,2, 3, 4. The paths with a single jump are given by IIo(r, . . . , r, r' , . . . r'), where 
(r, r') = (3,4) or (r, r') = (4,3), and where the jump takes place at location k = 
1, ... ,7V — 1. Note that the only jumps allowed are between r = 3 and r = 4, since 
T J (r,r') = if r ^ r' and r,r' {3,4} (see Proposition 3.2). 

It is thus clear that we have exactly 2(^) paths with k jumps. The factor 2 takes 
into account that r\ can take either of the values 3 or 4. For k fixed, the summand in 
(44) is bounded by 



| e iEf =1 ^-) (BVd.noCn, . . .,r N )AAc{) < C C 1 fc C 2 7V - fc , (108) 



where 



C = \\B*iPo\\ PH, 

C\ = maxmax|^(r,r')| =0(r max t/iV), 

C 2 = max max | T j (r, r ) | = 1 + O (r max t /N) . 
j r 

The last estimates on C\ and C2 follow from (107) and (106), respectively. The sum 
over all paths (N fixed) in (44) has the upper bound 

N-i //vN 

2C oE <^ k )C^- k = 2C [(C 1 + C 2 ) N -C( f ] <2C (l + 2CrUt/N) N , 

where C is such that C\ < Cr max t/iV and C 2 < 1 + Cr^ iax t/N. The limit as N -> 00 
of the r.h.s. 2Coe 2C,Tmax *. This implies that we can truncate in a controlled way the 
series over the number of jumps obtained in the continuous time limit. If we truncate 
at k < K, then the remaining tail of the series is estimated from above by 

N 

2C ( k ) C i C 2~ k <2C (C 1 /C 2 ) K+1 e 2CT ^ t = o((r^ x t) K+1 e 2C ^) . 

k=K+l ^ ' 

For bounded t, the tail of the series is thus 0((T max ) is ' +1 ) . 
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Contribution of paths without jumps. The products of transition coefficients for 
the constant paths with rj = 1,2 are 1. Also, 



N-l 



T 1 (3, 3) • • • r JV-1 (3, 3) = ] [ 



1 + 



y+Ut/N) 

l + y + (jt/N)' 



■(y + ((j + l)t/N)-y + (jt/N)) 



(109) 

and for r = 4 the product is given by the r.h.s. of (109) with y + replaced by y_. By 
taking the logarithm of (109) the product transforms into a Riemann sum, and so one 
easily obtains 



hm T 1 (3,3)-..r^ 1 (3,3)=exp{ f ds) 

U l + y+(s) 2 J 



! i + y + (t) 2 
i + y+(o) 2 ' 



The limit of the products with r = 4 is given by the latter square root with y + replaced 
by y_ . Finally, the limits of the exponential factors in (44) are 



lim e 1 ^ 1 '^^ = e'^o e(s ' rj)ds , 

N— >oo 



(110) 



where e(s,r) is defined as in (102). This gives the first line in (26). 

Contribution of paths with one jump. Let jo £ {1, . . . , N — 1} be the location 
of the jump. We have r\ = ■ ■ ■ = rj = r and r JO+ i = • • • = r^-i = r', where 
either (r, r') = (3,4) or (r, r') = (4,3). We treat the first the transition 3^4. The 
contribution to the sum in (44) is given by 



(B*Vo,Jl|??(0,3))(^,4)|^Vref), 



(111) 



where 



N-l 



30=1 



i^Ef =1 ^(3)+i^Ef = , 0+1 ^'(4) 



30-1 

n^(3,3) 

3=1 



T jo (3,4) 



N-l 

n ^(4,4) 

3=30+1 



(112) 



The continuous time limit of J\ is 

ft 



Jo 



y + (s)y'_(s) 
l+V+(s) 2 



w(s,t,4), 



(113) 



with w defined in (104). The corresponding quantity for the transition 4 — > 3 is obtained 
from (113) by interchanging the indices 3^4, replacing y + by y_ and y'_ by y' + . This 
completes the proof of Theorem 3.4. ■ 



3.2 Regime of isolated resonances 

In the regime of isolated resonances we can use Theorem 2.1. All that remains is to 
calculate the eigenvalues <P and eigenvectors rf and rp of the level shift operator (74). 
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Theorem 3.5 (Resonances) Suppose all gaps A 3 are numbers, well separated from 
0, independent of X. The resonances are 

e j (l) = , e 3 {2) = \)5 3 {2) + O(Aj) , (114) 

e*(3) = A 3 + \]S 3 (3) + 0{\j) , e 3 {4) = -A 3 + A 2 <P(4) + O(Aj) , (115) 

where, 

S j (2) = mj(A 3 ), (116) 



^) = ^ + 5^(^ I ),7(rAi)) I (118) 

7W = ^P / d CT | 5 (^H,a)| 2 (119) 
is 2 



2 

„2 



P (^t) ,7(rA^)^ = lim£ ^ X{ ~f}[ ^ 7 (rA^)dr (120) 

and i/ie eigenvectors are 

rf(l) = V++ + ^ , = — j ±—[ e -f>*<p ++ + ^__] , (121) 

rf(2) = ^ ++ - , 7j 3 (2) = -^L_ [ y++ _ e -^V__] , (122) 

^'(3) = »j?'(3) = , ^(4) = ^"(4) = . (123) 

The transition coefficients are T 3 (r,r') = (rp (r), rf +1 (r')), which equal 

T 3 \l,l) = T 3 (2,2) =T j (3,3) =T j (4,4) =1, (124) 

THl 2) = sinh ^ AJ+1 - AJ ]/ 2 ) (125) 
1 ' ; 2cosh(/3Ai/2)cosh(/3Ai+V2)' 1 ' 

All other transition coefficients vanish. 

Remarks. 1. If we consider the asymptotic regime where A 3 <C 1, we obtain 
e j (l) = 0, e 3 {2) = 2ia, e 3 {3) = A 3 + ia, e j {4) = -A 3 + icr. (To see this, note that 
(^-p 7(rA J )) = 0((A 3 ) 2 ) in this regime.) This agrees with Theorem 3.1 if one takes 
t 3 ' ^ oo, which signifies A 2 < A' < 1. The biggest change is the change to t 3 , which 
previously was just A 3 /a 3 . Now we have 

f 00 1 - X(-ee)(r 2 - 1) 
a(A 3 )r(A 3 ) = A 3 - lim J J_\ L l{rA 3 )dr. (126) 

The principal value integral appearing in (126) vanishes in the limit A 3 — > 0, and is 
therefore part of the remainder in the setting of Theorem 3.1. 
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2. Clearly the full analysis of the overlapping region is more involved than taking the 
limit of the answers from the non-overlapping region. For instance, the limit above does 
not give the correct answer for the eigenvectors, and hence nalso for the transmission 
coefficients. 

3. Once Theorem 3.5 is proved, the calculations leading to (29) and (30) for the 
single sudden crossing are done exactly in analogy to Section 3.1.1, and in fact are 
easier. 

Proof. We have to calculate, and diagonalize the level shift operator (74). For this 
purpose we will take j to be fixed, in order to prove (114) - (123). (The calculation 
of the transition coefficients follows trivially from these, using the definition (87).) We 
will consider the gap A J to be positive. Then I = V — </rA r ' VA r Jr. Using (147), 
(148) and (150), we see that 



(r p g)(-l) N + (-ira{e-^T fi g) 



I = a x <g> Is <8> <t>{rpg) - l s <8> cr x & 
and the level shift operator at energy e is 

A(e) = P e IPe(L R + L s -e + iO)- 1 P e IP t 



N , 



(127) 



(128) 



with L s = — (cr z <8> Is - Is ® a z) ® 1r and Lr = Is <8> Is <8> dr(u). Also P is the 



projection onto the span of {(f++ <&£l,(p ® f2}. So, in this basis 

A(0) 
where 

2A n (0) = (rpg, (u - A + iO)" 1 ^) + (e~^g, (u + A + iO)" 1 ^ 



'An(O) A 12 (0) 
A 2 i(0) A 22 (0)_ 



2A 22 (0) = (rpg, (u + A + iO)" 1 ^) + (e 



-(iu 



A + iO)- 1 ^ 



-2A 21 (0) = (e-^rpg, (n - A + iO)" 1 ^) + ( Tp g, {u + A + iO)" 1 ^) 
-2A 12 (0) = (e-P u Tpg, (u + A + iO)" 1 ^) + ( Tp g, (u - A + iO)" 1 ^) 



where the inner products are in L 2 (R x S^dudo-). Noting that e-^rpg^a 
\[Tpg](— u,a)\ 2 , we see that 

An(0) = -A12 = i<T /J 5,Im[(«- A + iO)- 1 ]^) , 
A 22 (0) = -A 21 = i ( Tp g, Im[(n + A + iO)" 1 ]^) . 



(129) 

(130) 
(131) 
(132) 
(133) 

| 2 = 

(134) 
(135) 



But now we use the well-known formula that lim e j Im[(u — uq + ie) l ] = — it5(u — no), 
as a distribution to obtain 



A(0) = -in [ 
Js 2 



\r p g(A,a)\ 2 -| W (A,a)| 2 ' 
-| W (-A,a)| 2 \rfsg(-A,a)\\ 



da 



i7T7(A) 
e-^ + l 



-1 



(136) 
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It is easy to see that the eigenvalues are 5(1) = and 5(2) from (116), with eigenvectors 
r/(l), r/(2), r/(l) and 7/(2) from (121) and (122). 

For e = ±A there is only one eigenvector, each. Let us consider e = A whose 
eigenvector (left and right) is ip + - . In this case the action of a x <8> Is and Is ® 0":r both 

serve to map to eigenvectors of Ls with eigenvalue 0. On the other hand, e = A 

appears in the resolvent. So 

2A(A) = ( Tp g, (u - A + iO)" 1 ^} + (e^Tpg, (u - A + iO)" 1 ^) 
= (r^, [(« - A + iO)- 1 + (-u - A + iO)" 1 ] r^> . 

(We again used e _/3u | [rg<7](u, a)\ 2 = |[t/3p](— u,ct)| 2 .) Now we use the well-known for- 
mula that lim e ^o(u— uo + ze)" 1 = — i7rJ(n — tto ) + V ( u } ) , where V( u ^ ) is the Cauchy 
principle value distribution, V( u } uo ) = lim e _>o X| ^ t l^ >e , where the limit is in the space 
of distributions. So in this case we do get a "Lamb shift" in addition to the purely 
imaginary resonance 

2A(A) = { Tp g, [-m[5(u - A) + 5(u + A)] + V(u - A) - V(u + A)]r^g) . (138) 

Consideration of this formula leads to (117), and (118) then follows by symmetry 
arguments. 

We can treat the case of negative gaps A" by conjugating by a x to change to — A J . 
This does not affect the resonances but it does affect the eigenvectors. However, all that 
happens is that some vectors rj J (r) are multiplied by —1 and this is always accompanied 
by the same change to the corresponding dual eigenvector 7p(r). This type of gauge 
transformation does not affect any physical quantities. □ 



A Araki-Wyss representation 
A.l The representation 

Here we will outline the Araki-Wyss representation [5] in order to be self-contained. 
We use similar notation to [10] (starting on page 24). The Araki-Wyss representation 
of the CAR is a representation on the tensor product of two fermionic Fock spaces 
JF(L 2 (M 3 , d 3 k)) <8> JF(L 2 (M 3 , d 3 k)), such that the smoothed-out creation operators are 
represented by the formula 

= ( v^TT 9 ) 1 + M) " ° ( v^TT 5 ) ' (139> 

where N is the number operator A" = dT(l). The annihilation operators are ap(g) = 
[a*p(g)]*. Let us write (ip = (1 + e _/3 ' fc ' In this context, the modular operator and 
modular conjugation combine to give 

ap(g) = JnA^ 2 ap(g)A^ 1/2 J R 

( e^ fc l 2 /2 \ „ M „ / e /3|fc| 2 /2 \ (140) 
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and 

a}{g) = J R AK 2 al(g)A R 1/2 J R 

( e -/3|fe| 2 /2 \ / e -/3|fc| 2 /2 \ M (141) 

VVe^+1 / \\/rW + 1 / 

A. 2 Regularity of form factors 

Let (/ £ L 2 (R + x S" 2 , |/c| 2 d|/c| da) be a form factor represented in spherical coordinates 
(da being the uniform measure on S 2 ). We introduce a new radial coordinate r = \k\ 2 
so that the dispersion relation of the Fermions reads u> = \k\ 2 = r, i.e., H R = dT(r) 
on the Fock space J r (L 2 (R + x S 2 , -^drdcr)) (see (4)). The Araki-Wyss representation 
Hilbert space associated to the thermal equilibrium is 

F(L 2 (R + x S 2 ,^drda))®F(L 2 {R + x S 2 ,^drda)). (142) 
For the purpose of spectral deformation, it is advantageous to use the maps 

fl # (/)®lwa # (/9 0) (-1) N ®a*{f) ^ a # (0© /) (143) 

to define an isometric isomorphism between (142) and the Hilbert space 

f(l 2 (R + x S 2 , ^dr da) © L 2 {R + x S 2 , ^dr d<r)) . (144) 

In (143), N is the number operator. A further isometric isomorphism between (144) 
and 

H R = f(l 2 {Rx S 2 ,duda)) (145) 



(see also (32)) is induced by such an isomorphism between the one-particle spaces, 
given by 

i- ^ 7 | / \ M 1 / 4 f f(u,a) if u > 0, M ,.s 

V ; v/2 I 0(-«,<7) if «<0, v ; 

where f,g £ L 2 (R + x S^^drda) and /i G L 2 (M x S 2 ,dudcr). 

Under these isomorphisms, the field operator 4>(g), in the Araki-Wyss representa- 
tion, has the expression 

Md) = tlrpg), (147) 
where <fi on the r.h.s. is the field operator on the fermionic Fock space (144), and where 



M(^-) = ^y^ TT {^/; )) ifn - ; d48) 

Using the explicit transformations introduced above, the expressions (37)- (41) are easily 
found. Furthermore, it is not hard to analyze the explicit action of the deformation 
transformation Ug (see before (50)). For instance, for h € L 2 (R x S 2 ,duda) and 9 6 M, 
we have 

Ug<t>{h)U* e = <f>(e dd "h), (149) 
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where (e S9u h)(u, a) = h(u+6, a). One can now combine (149) with h = Tpgi (see (147), 
(148)) and check that due to Condition (R) in Section 1 (with a(u) = 1), 6 i-> Ue4>(h)Ug 
admits an analytic continuation into a strip |Im#| < 9 . One can proceed similarly to 
check analyticity of (V^){6) by using the explicit form Ar = e _/3LR for the modular 
operator of R. For reference in the body of the paper, let us note explicitly that 

J R Ajf <^(<?) J R Ajf = a*{T P g)(-\) N + {-\) N a^rpg) , (150) 

in this representation .F(L 2 (]R x S 2 ,duda)). This can also be seen by comparing to 
(140) and (141). 
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